Abundant crystalline silicates in the disk of a very low mass star 
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§ ■ ABSTRACT 

We announce the discovery of SST-Lup3-1, a very low mass star close to the brown dwarf 
boundary in Lupus III with a circum(sub)stellar disk, discovered by the 'Cores to Disks' Spitzer 
Legacy Program from mid-, near- infrared and optical data, with very conspicuous crystalline 
| silicate features in its spectrum. The Spectral Energy Distribution of the system is fit with 

J> ■ stellar and disk models to study the disk structure. SST-Lup3-1 is a M5.5 star with a luminosity 

t-H I of 0.08 L Q (at an assumed distance of 200 pc), a mass of 0.10 Mq, a radius of 1.09 Rq, and an 

age of ~ 1 Myrs. It is the first of such objects with a full 5 to 35 fim spectrum taken with the 
IRS and it shows strong 10 and 20 fim silicate features with high feature to continuum ratios and 
clear crystalline features out to 33 /Ltm. The mass of its flared disk is constrained to lie between 

[ — | 2xl0 -4 and 10~ 7 M Q (80-0.03 M ffi masses) and the dust in the upper layer has a crystalline 

silicate grain fraction between 15% and 33%, depending on the assumed dust continuum. The 
availability of the full Spitzer infrared spectrum allows an analysis of the dust composition as 
a function of temperature and position in the disk. The hot (^300 K) dust responsible for the 
10/im feature consists of a roughly equal mix of small (~ 0.1 /zm) and large (~ 1.5 fim) grains, 
whereas the cold (~70 K) dust responsible for the longer wavelength silicate features contains 
\ primarily large grains (> 1/im). Since the cold dust emission arises from deeper layers in the 

inner (< 3 AU) disk as well as from the surface layers of the outer (3-5 AU) disk, this provides 
direct evidence for combined grain growth and settling in the disk. The inferred crystalline mass 
fractions in the two components are comparable. Since only the inner 0.02 AU of the disk is warm 
enough to anneal the amorphous silicate grains, even the lowest fraction of 15% of crystalline 

■ material requires either very efficient mixing or other formation mechanisms. 

Subject headings: circumstellar disks - very low mass stars - stars: individual (SST-Lup3-1) 
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1. Introduction 

Ground based observations of young brown 
dwarfs in star forming clouds have shown the 
presence of significan t infrared (IR) excesses (e.g. , 
Comeron et al. 2000l ; Javawardhana et al. 2003 ) 



and broad Ha emission l ines (e.g.. lMuzerolle et al 



2003; iNatta et alJ 12004 ) suggestive of the pres- 
ence of circum(sub)stellar disks and accretion in 
these very low-mass objects. Also, recent ground- 
based mid-IR photometry has shown hints of grain 
growth and dust settling from the silicate emis- 



sion f eatures of some of these disks (jSterzik et al 
2004t lApai et all l2004h . Observations with the 



Spitzer Space Telescope (SST or Spitzer) permit 
the first d etailed studies of the geometry of such 
disks ( e.g., Allers et alJl2009) and of their silicates 
(e.g., iFurlan et al.ll2005l )~ which provide interest- 
ing comparisons with the dust evolution in disks 
around higher-mass T Tauri and Herbig Ae/Be 
stars. 

It is known that the presence and intensity 
of amorphous silicate features in T Tauri stars 
anti-correlates with the size of the gra i ns in 
the disk surface layer s (iBouwman et al. 2001 
van Boekel et all 120031 . 12005^ IPrzygodda et al 



20031 iKessler-Silacci et al.1 120091 Upai et al 



(|2005l ) analyzed Spitzer spectra of the 10 /im sili- 
cate feature of 6 brown dwarfs and very low mass 
stars with disks in Chamaeleon I and found an an- 
ticorrelation between the crystallinity of the sili- 
cates and the degree of disk flaring, which suggests 
an evolutionary process in which the disk flattens 
and the dust crystallizes as the grains grow and 
settle to the disk mid-plane. This explanation, 
although very attractive, is based on observations 
of only the 10 /im silicate feature, with low peak 
to continuum ratios and with spectra dominated 
by either crystalline or amorphous features but 
not in combination. 

In the present paper we announce the discov- 
ery of SST-Lup3-1, a new very low mass star close 
to the brown dwarf boundary with a disk in the 
Lupus III dark cloud, and analyze the properties 
of its disk from its Spectral Energy Distribution 
(SED). Its Spitzer InfraRed Spectrograph (IRS) 
spectrum over the entire 5-35 /im range shows 
strong amorphous and crystalline silicate emission 



features and is presented and analyzed here to ex- 
tract the mineralogical composition of the dust in 
the disk. 

2. Target selection and observations 

Th e 'Cores to Disks' (c2d) Spitzer Legacy pro- 
gram (jEvans et al.l 120031 ) has mapped five nearby 
large molecular clouds with IRAC and MIPS at 
6 wavelengths between 3.6 and 70 /im. Photome- 
try was extracted for all point sources and band- 
merged with the 2MASS catalog. For more in- 
formation s ee the Delivery Do cumentation of the 
c2d project ( Evans et al. 2006h . The area has also 
been surveyed with the 2.2m telescope of the Euro- 
pean Southern Observatory (ESO) and the Wide 
Field Imager (WFI) instrument in the R, I, and z 
broad-band filters to complement the Spitzer data 
and study the very low-mass stellar population 
(Comeron et al. in prep.). Detailed descriptions 
of the observations of Lupus by the c2d Team can 
be found in Mundy et al. (in prep.) for IRAC and 
in lChapman et all (|200d ) for MIPS. 



Using the band-merged catalog of the Lupus 
III dark cloud, we have discovered a rich popula- 
tion of IR excess sources with luminosities in the 



substellar domain according to the iBaraffe et al 
(|1998l ) pre-main sequence tracks. Constructing 



NTT at La Silla (Chile) for program 276.C5031. 



their SEDs with WFI RIz, 2MASS plus IRAC 
and MIPS fluxes from 0.6 to 70 /im and com- 
paring them with stellar emission models, we find 
that some of the IR excesses at near and mid IR 
wavelengths may be explained with the presence 
of circum(sub)stellar disks. We selected one of 
them with optical and near IR colors characteris- 
tic of a late M-type dwarf which could be observed 
with the Spitzer IRS. The object, with coordi- 
nates 16 h ll m 59.9 s , -38°23'37.5" (J2000), is named 
SSTc2d J1 61159.9-382337. We na me it SST-Lup3- 
1 following [Comeron et all (|2003h . and nickname 
it Veronica's star. 

SST-Lup3-1 was observed with the IRS on- 
board Spitzer on 14-08-2005, with AOR Key 
0015737856, in both short low (SL) and long low 
(LL) modules. The total wavelength coverage is 5 
- 35 /im and the integration times were 120 and 
620 s for SL and LL, respectively. Data reduction 
started from the Basic Calibrated Data images, 
pipeline version S12.4.0. The processing includes 
bad-pixel correction, extraction, defringing, and 
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order matching using the c2d analysis pip eline 



(|Kessler-Silacci et al.ll20o5 lLahuis et al.ll2006h 



In addition, an optical spectrum of SST-Lup3-1 
was obtained with the EMMI instrument on ESO's 
New Technology Telescope. The object was ob- 
served with the grism # 6, which provides a spec- 
tral resolving power of 1,500 with the 1" slit, and 
with 9 minutes of integration time. The spectrum 
was reduced using standard IRAF tasks for extrac- 
tion, wavelength calibration, and relative flux cal- 
ibration within the 'onedspec' package. Relative 
flux calibration was performed by using the spec- 
trum of the spectrophotometric standard 8 Vir as 
a reference. 

3. Stellar parameters 

Figure 1 shows the full optical spectrum of 
SST-Lup3-1 compared to spectra of young M-type 
stars. It is clear that the object belongs to an in- 
termediate class between the M5 and M6 spectral 
types and that the moderate extinction towards 
the source makes it resemble a colder M6.5 brown 
dwarf. The TiO band-head absorptions at 6200, 
7100 and 7600 A (|Kirkpatrick et al.lll99lh and the 
slope of the continuum around 7500 A were used 
for the spectral typing after the extinction correc- 
tion. 

We subsequ ently used an impro ved version of 
the method by iFrasca et alJ (|2003l ). in which the 
spectrum is compared with a grid of standard 
spectr a of M to L dwarfs from iKirkpatrick et al" " 



(|199ll ) with different extinctions to get the best fit 
to both parameters simultaneously. This yields a 
spectral type of M5.5 with an error of less than 0.5 
spectral subclasses and an extinction of Ay=2.2 
mag. A similar analysis with the method by 



Guieu et al 



(2006), which uses instead an aver- 
age of dwarf and giant M to L field standards to 
better approach the emission of young stars, yields 
a spectral type of M5.75 and Ay— 2.08 mag. The 
similarity of these two results gives confidence in 
our determination of the spectral type of the ob- 
ject. 

The spectrum of SST-Lup3-1 shows Ha (6563 
A) in emission with an equivalent width of 16.2 
A. The profile is too narrow (FWHM ~ 5 A) to 
be a signature of disk accretion an d is most likely 
related to chromospheric activity (jWhite fe Basri 
2003, Ijavawardhana et al.l 120021 ). However, its 



presence in the spectrum, as well as in the other 
M6 and M6.5 young stars shown in the figure, in- 
dicates the youth of the object and confirms its 
membership to the Lupus molecular cloud. The 
detection of the Li I absorption line at 6708 A with 
an equivalent width of 0.5 A supports the young 
nature of the object. 

The SED of the object is constructed using the 
photometry in Table 1 and compared w ith Star- 
Du sty stellar models dXllard et alj |2000h. Follow- 
ing Luhman et al. ( 20031 ). a young M5.5 star has 
an effective temperature of 3057 K and intrinsic 
color R c — I c = 1.88. Given the observed R c — I c 
color of 2.00 we derive an extinction from the opti- 
cal ph otometry of A v = - 67 us ing the extinction 
law of iRieke fe Lebofskvl (|1985f ) . The difference 
in the extinction determinations from the optical 
photometry and spectrum is larger than their un- 
certainties and cannot be due to the object hav- 
ing surface gravity of a giant or dwarf. However 
it would disspaear if the star would have an M5 
spectral type, in which case A v ~ 2, consistent 
with the spectroscopic result. Possible reasons for 
this discrepancy are the scatter in the intrinsic 
R c — I c colors of M-type stars, that the object 
is variable at optical wavelengths, as seen in other 
very low mass stars (jCaballero et al.ll2004l ) , that it 
presents bluer optical co lors than those expec ted 
for its spectral type as in lLuhman et al.l (1998), or 
that the standard interstellar extinction law is not 
valid for the Lupus dark cloud. 

The difference in extinctions introduces only 
small differences in luminosity of 0.015 L Q and 
age of 0.4 Myrs so it is not critical for the subse- 
quent analysis. We thus adopt an extinction value 
of A v — 0.7, which is consistent with all of the 
optical and near-infrared photometry, and gives 
a luminosity of 0.081 L Q , or log L/L = —1.09, 
for an assumed distance to Lupus III of 200 pc 
(|Comer6rJl2006h . Together with the T off of 3050 
K, this gives an age of ~ 1 Myr, a m ass of 0.10 
and a radius of 1.09 Rq using the iBaraffe et al 



(1998) isochrones. 



4. Disk parameters 

The aim of this work is to link the analysis of 
the SED, which is related to the geometry of the 
disk, and the mid-infrared spectrum, which gives 
the mineralogical composition of the dust in the 
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6000 



Fig. 1 — NTT/EMMI optical spectrum of SST- 
Lup3-1 compared to spectra of M4 (Gl 402) 
and M5 (GJ 1028 ) spec troscopic standards from 
Kirkpatrick et al.l (|l99lh and M6 (Sz 112) and 
M6.5 (Sz 109) young dwarfs in Lupus III from 



Comeron et al 



(12003). The thin line corresponds 
to a dereddened spectrum with Ay— 2.2 and the 
dashed line is the observed spectrum. The upper 
inset shows the Li I absorption line at 6708 A. 



Table 1: Photometry of SST-Lup3-1 



Band 


A 


Magnitude 


Flux 


Rcf. 




/im 


mag 


mjy 




R 


0.70 


16.00 ± 0.10 




1 


R, 


0.64 


16.251 ± 0.009 




2 


Ic 


0.79 


14.256 ± 0.007 




2 


z 


0.96 


13.597 ± 0.008 




2 


J2MASS 


1.25 


12.197 ± 0.024 




3 


H2MASS 


1.62 


11.511 ± 0.026 




3 


K2MASS 


2.20 


11.204 ± 0.023 




3 


IRAC-1 


3.60 




17.03 ± 2.55 


4 


IRAC-2 


4.50 




14.60 ± 2.19 


4 


IRAC-3 


5.80 




11.70 ± 1.76 


4 


IRAC-4 


8.00 




12.90 ± 1.94 


4 


MIPS-1 


24.0 




22.40 ± 3.36 


5 


MIPS-2 


70.0 




< 50.00 
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Fig. 2.— Spectral Energy Distribution of SST- 
Lup3-1, discovered in Lupus III with the c2d 
Spitzer data. The photometry is described in Ta- 
ble 1. Open and filled circles are the observed 
and dereddened fluxes respectively. The dotted 
line is the stellar model of a M5.5 dwarf. The full 
Spitzcr/IRS spectrum of the object is overplotted 
on the photometry and the total flux is fitted with 
a low mass CGPLUS flared disk model called 'high 
continuum' (dashed line) and a heavier D'Alessio 
et al. flared disk model called 'low continuum' 
(dot-dashed line) (see text). The inset illustrates 
the three different continua used to analyze the 
IRS spectrum; the full line is the usual power-law 
continuum. 



References: 1. DENIS catalog (2005), 2. Comeron 
et al. (in prep.), 3. Cutri et al. (2003), 4. Mundy 
et al. (in prep.), 5. Chapman et al. (2006), 6. 
Young et al. (2005) 
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upper disk layer. The problem has no unique solu- 
tion because the dust composition is commonly de- 
termined by subtracting the disk continuum SED 
from the IRS spectrum whereas the disk geometry 
and the SED itself are determined partially by the 
dust properties. A larger sample of ~ 100 spec- 
tra of disk sources will be analyzed with the same 
method (Olofsson, Augereau et al. in prep.); here 
we introduce the method and the relevant ranges 
of disk and dust parameters for this specific object. 

We have computed flared irr adiated disk mod- 
els with the CGPLUS code ( Dullemond et al 



200 ll ) and flared irradiated accretion disk models 
with the prescriptions of D'Alessio et al. (|2001 



2005) using the stellar parameters of SST-Lup3-1 
as inputs for the central star. The dust compo- 
sitions were chosen to produce negligible silicate 
emission: pure carbon grains of interstellar size 
in CGPLUS and large grains in the disk surface 
in the D'Alessio et al. models. This provides 
possible disk continua to subtract from the IRS 
spectrum assuming that the silicate features in the 
spectrum are produced in an optically thin disk 
upper layer. A model is considered a good fit for 
this analysis when it fits the dereddened 2MASS 
and IRAC photometry, matches the IRS spectrum 
shortward of 8 fim and does not exceed it at longer 
wavelengths. This yield a range of disk parame- 
ters that fulfill that condition: the CGPLUS disk 
models give an upper limit to the disk mass (gas 
plus dust with the usual ratio of 100:1) of 10~ 7 
M , inclination angle i > 70° (90° being edge-on) 
and i?disk < 15 AU. The D'Alessio et al. models 
provide reasonable fits for disk masses between 
2xl0~ 4 and 0.018 M Q , inclinations i > 30° and 
^disk < 50 AU. The large difference in disk prop- 
erties in the computed cases comes mostly from 
the different dust opacities in the two sets of mod- 
els: the D'Alessio et al. models allow a range of 
maximum grain sizes from 1 /im to 10 cm, which 
make the disk opacities low enough to get very 
small infrared fluxes with larger disks than the 
CGPLUS models, which use the dust opacity of 
astrono mical carbon and silica te sub- micron sized 
grains (jDullemond et al.l 120011 ). In any case, the 



quoted disk radii and masses of these disk models 
are representative only of the material probed by 
the IRS spectrum and therefore lower limits to the 
actual values. Both sets of models produce flared 
disks with H(R)/R oc R 2 and E(R) oc R' 1 . 



Figure 2 shows the SED of SST-Lup3-1 to- 
gether with example models, with properties con- 
veniently chosen to explore the range of reason- 
able continua for the subsequent analysis of the 
IRS spectrum. The dot-dashed line corresponds 
the 'low continuum' model, a flared D'Alessio et 
al. disk model with an outer disk radius of 10 AU, 
an inner radius of 0.01 AU (3 times the stellar ra- 
dius), a mass of 2.4 xl0~ 4 M Q , a mass accretion 
rate of 10~ 9 M Q yr _1 , a viscosity a = 0.01, a max- 
imum grain size of 100 fim and an inclination of 
60°. It represents a relatively massive disk with 
low dust opacities in this comparison. The dashed 
line is the 'high continuum' model, an extremely 
low mass CGPLUS disk model with an outer ra- 
dius of 15 AU, an inner radius of 0.06 AU (20 
times the stellar radius), a mass of 10 -7 M and 
an inclination angle of 75°. More massive disks 
with these small dust grains rapidly produce larger 
fluxes than those seen in the IRS spectrum. The 
inset shows a blow-up of the IRS spectrum with 
both continua together with the usual power-law 



fit to the continuum (as in Kessl er-Silacci et al 



2006) for comparison. Both disk models (very low 
mass with very small dust particles or more mas- 
sive but with large dust particles) produce similar 
SEDs up to the mid-IR wavelengths but have very 
different shapes at long wavelengths. The need for 
large grains or very small masses of small grains 
suggests that some degree of dust grain growth 
and settling to the midplane have taken place in 
the disk of SST-Lup3-1. 

5. Dust composition in the disk 

Figure [3] shows the IRS spectrum of SST-Lup3- 
1 from which the disk models in Figure 2 have been 
subtracted. Both spectra are then compared with 
the best-fit silicate emission model. The spectrum 
has, to our knowledge, one of the highest peak to 
continuum ratios (> 2.0) in the 10 /im silicate fea- 
ture ever published for a disk around a very low 
mass star. The spectrum also shows evidence of a 
moderately high degree of crystallization, with a 
Sn.s/Sg.s ratio equal to 0.9, and conspicuous fea- 
tures from crystalline silicates such as forsterite a t 
11.4, 23.8, 27.9, and 33.7 fim ([Koike et al.lll993h . 
The 1-<T error bars in the upper panel quantify 
the signal to noise around these features. Prob- 
lems with the extraction of the third order of the 
LL module at wavelengths between 19.5 and 21.5 
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/zm for very faint sources do not allow the identi- 
fication of the 19.8 or 21.0 /im forsterite features. 
The 10.2 and 16.5 /im features are not significant 
within the noise. The rest of the crystalline fea- 
tures are present in both full ape r ture a nd profile 
fit extractions, see Lahuis et al.l ( 20061 ) for more 
information about the quality assessment of the 
c2d IRS data reduction pipeline. 

5.1. Method and results 

Compositional fits to the full IRS spectrum 
of SST-Lup3-1 have been obtained fo l lowin g: the 
method outlined in Ivan Boekel et a l. (2005) and 



used in lApai et al.l (|2005l ) and in iHonda et al 
( 20061 ): a x 2 minimization between the spectrum 



and a weighted sum of dust opacities, multiplied 
by one or more blackbodies, for the five principal 
dust species and for 2 grain sizes (0.1 /im and 1.5 
yitm). For every fixed blackbody temperature, this 
procedure provides relative mass abundances, and 
a minimum \ 2 value. This results in a relation 
between the temperatures and the minimum x 2 
values, from which we infer the best-fit temper- 
ature. We use the distribution of hollow spheres 
routine for ca lculating the op acities of the crys- 
talline grains ( Min et al. 20051 ) and the standard 
Mie theory for the amorphous grains. 

Two novelties are applied here: the first one 
is that we use disk models to estimate the disk 
continuum emission, which connects the composi- 
tional analysis with the physical picture of the disk 
and allows us to study the height and/or depth in 
the disk of the dust emitting the 10 and 20 /im 
amorphous and longer wavelength crystalline sil- 
icate features. The second one is that we fit the 
entire IRS spectrum from 5-35 /im, which re- 
quires the inclusion of an extra low-temperature 
dust component for the long wavelength part and 
thus allows different parts of the disk to be probed. 
This long wavelength range is also essential to as- 
certain the presence of crystalline silicate features, 
many of which arc potentially confused with PAH 
features shortward of 13 /im. 

We first fit the 10 /zm silicate feature consid- 
ering only the spectral range between 7.5 and 
13.5 /im, which is least affected by the continuum 
choice and which requires a relatively high temper- 
ature of ^300 K. This best fit is then subtracted 
from the spectrum and the residual spectrum is fit 
with a lower temperature component with a best- 
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Fig. 3. — Continuum-subtracted Spitzer IRS spec- 
trum of SST-Lup3-1 and compositional analysis 
of the silicate emission. The thin line is the IRS 
spectrum of the source from which the low and 
high continuum disk models have been subtracted 
in the upper and lower panels respectively. The 
tickmarks identify the forsterite emission features. 
The upper panel shows the 1-cr error bars. The 
best-fit models are shown in dot-dashed lines and 
combine the emission of amorphous and crystalline 
silicate grains with sizes of 0.1 and 1.5 /zm of two 
different temperatures quoted in the plot. The fits 
show larger crystalline fractions for the case with 
high continuum and imply a minimum of 15% of 
crystalline grains in the disk. 
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fit value of ~70 K. This method has the virtue 
that it allows a direct comparison with fits per- 
formed to only the 10 /im feature (as in Apai et 
al. 2005, van Boekel et al. 2005), where only the 
high temperature component is considered. It also 
allows a direct determination of the mass fractions 
in each component, which are 59 and 41% for the 
hot and cold components respectively if the low 
continuum model is used. With the high contin- 
uum, these numbers change to 65 and 35%. Also, 
the two temperature approach identifies (via the 
X 2 minimization) the characteristic temperature 
ranges of the dust producing the amorphous 10 
and 20 /im silicate features. 

Table 2 shows the abundances of the different 
dust species that result in the best fit to the IRS 
spectrum shown in Figure 3 after subtracting the 
low continuum disk model. The abundances are 
found by averaging a number of realizations of the 
fit in which the spectral ranges used for the fits 
were changed and in which a range of tempera- 
tures was explored. The associated uncertainties 
are therefore larger than indicated by the stan- 
dard deviations of the data, since they account for 
the relatively low number of realizati ons. How- 
ever, t hey are similar to those cite d in lApai et al 



(|2005l ) and in Irlonda et all (|2006f ) in their fits to 
the 10 /im spectral range (high temperature com- 
ponent). The uncertainties in derived abundances 
are slightly larger for the longer wavelength spec- 
tral range. The quantitative fit results obviously 
depend sensitively on the adopted continuum, but 
the qualitative results are very similar. Hence, we 
will discuss primarily the low continuum case be- 
low, and mention the high continuum case only 
when it gives a different trend. Also, Table 2 
only provides a combined mass percentage for the 
amorphous silicates, which includes olivine, pyrox- 
ene and silica, because the fitting method was not 
sufficiently robust to identify the individual amor- 
phous species responsible for the long-wavelength 
part of the spectrum. This docs not affect the rel- 
ative mass fraction of amorphous and crystalline 
material, which was always consistently recovered. 

The best fit composition shows that i) the 
amorphous hot silicates are a roughly equal mix 
of small and large grains, which is also found in 
other sources with high peak to continuum ra- 
tios (~2.0 ) of the 10 Mm silic ate feature (see e.g. 
Fig. 3 in iHonda et al.1 [20061 ); ii) the crystalline 



silicates in the hot component are a mix of large 
and small grains independently of the chosen con- 
tin uum, a result also cons istent with that given 
by Ivan Boekel et~aL (2005) for similar 10 /im fea- 
tures; iii) the amorphous cold grains are mostly 
large (1.5 /im). This is a new and surprising result. 
However, this conclusion depends on the adopted 
continuum, since the high continuum yields mostly 
small grains in the cold component; iv) the cold 
crystalline grains are found to be smaller than 
2 /im regardless of the continuum used. 

In the following, we discuss in detail the total 
crystalline mass fraction and the distribution of 
large and small grains in the disk. 

5.2. Crystalline mass fraction 

The best-fit crystalline mass percentage varies 
from 33% for the high disk continuum, through 
20% with the power law continuum and down to 
15% with the low disk continuum. In general, 
it is found that the smaller the peak to contin- 
uum ratio, the larger the mass crystalline frac- 
tion. The case with the low continuum based on 
the D'Alessio et al. (2005) models represents the 
best fit to the spectrum, mostly at 10 /im, and also 
the lower limit to the total dust crystalline frac- 
tion, therefore all subsequent analyses will focus 
on this case. 

The percentages of crystalline grains in the hot 
and cold components are 16.7 and 12.2 % respec- 
tively, with a similar ratio but larger absolute val- 
ues in the high continuum case. The clear iden- 
tification of the 23.7, 27.9 and 33.7 /im forsterite 
features in the long wavelength spectrum is a good 
assessment of the presence of crystalline grains in 
the disk and confirms the result at short wave- 
lengths. We have a lso co mpared our results to 
those of Apai et al. ( 20051 ) for sources with sim- 
ilar 10 /im feature shapes using the high contin- 
uum case, which is the closest approximation to 
the definition of the continuum of those authors. 
Our inferred crystalline mass fractions are similar 
to better than 10%. 

5.3. Distribution of large and small grains 

Regardless of the adopted continuum, our fit- 
ting technique of the hot component recovers the 
approximately 50-50% of small (0.1 /zm) and large 
(1.5 /zm) grains found by previous analyses for this 
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Table 2: Mineralogical composition of the disk around SST-Lup3-1 after subtraction of the low continuum 
disk model. 



Species 


Mass % a 




Mass % a 






284±27 K 




73±9 K 






0.1 /im 


1.5 /im 


0.1 /im 


1.5 /im 


Amorphous silicates 6 


33.5±2.7 


49.8±5.3 


0.0+1.5 


87.8+9.6 


Crystalline forsterite 


6.0±0.5 


2.3±0.2 


11.8±1.2 


0.2+0.1 


Crystalline enstatite 


0.7±0.1 


7.7±0.6 


0.0+1.5 


0.2+0.1 



a Relative fractions of 59% (284 K) and 41% (73 K). 
b Includes olivine, pyroxene and silica. 



type of 10 /im feature. The fit to the long wave- 
length part of the spectrum yields mostly large 
grains for the low continuum (87.7 % of large 
grains, Table l5~Tj) and mostly small grains for the 
high continuum (76.3 % of small grains). 

To test the sensitivity of the results to the 
adopted large grain size, a series of fits to the long 
wavelength range of the IRS spectrum was per- 
formed with 1.0 and 5.0 /im grain sizes but this 
did not improve the fits dramatically. The re- 
sults were consistent with what has been reported 
above, namely that most of the grains in the 
low-temperature component are 1 /im-sized amor- 
phous grains for the low continuum case. There- 
fore, even with the small difference in feature pro- 
files for grains in the 1-5 /im size range, the best 
fit always selects grains of ~l/im radius for the 
low temperature component. 

The crystalline grains are also a mix of large 
and small grains in the hot component regard- 
less of the continuum used. However, in the tests 
with larger grain sizes for the cold component, we 
find that the best fit always selects the smallest 
available grains, with grains of 0.1 - 1.5 /im giv- 
ing the best fit to the long wavelength spectrum. 
In these cases, the cold temperature components 
of the silicate emission model (dashed lines in Fig. 
3) show the ~1.T:2 relative flux ratios between the 
23.7, 27.9 and 33 .7 /xm forsterite features fo r small 
grains (Fig. 8 in Kessler-Silacci et al. 20061 ). 



6. Discussion 

In this section we connect the results of the 
compositional analysis with the disk modeling to 
provide a combined interpretation of all the data. 
As shown in Table 2, the hot temperature emission 
component contributes 59% of the total dust mass, 



reproduces completely the 10 /im feature and con- 
tributes to the 20 /im feature. The cold tempera- 
ture component, mostly coming from large amor- 
phous grains, reproduces the rest of the 20 /im fea- 
ture and accounts for most of the forsterite fea- 
tures longward than 20 /im. Figure [?] shows the 
2D temperature structure of the D'Alessio et al. 
disk model (from Sect. 4) used to fit the entire 
SED except the silicate emission features. The 
~300 and ~70 K contours indicate the locations 
in the disk surface where the different features are 
mostly produced. The hot inner disk region (from 
0.05 to 1 AU) with temperatures around 300 K 
emits mostly in the 10 /im range. The longer wave- 
length emission comes from dust at a temperature 
of around 70 K, which is found both in the outer 
disk surface (from ~3 to 5 AU) and in the deeper 
disk layers at smaller radii. 

Interestingly, the majority of the grains in the 
low temperature component are large grains, while 
only half of the grains in the high temperature 
component are large. Together with the distribu- 
tion of dust at 70 K in the disk, this strengthens 
the possibility that part of the long wavelength 
emission spectrum comes from deeper layers of 
the inner disk rather than just the outer disk sur- 
face. This could be interpreted as a signature of 
sedimentation to the disk mid-plane. This con- 
clusion is not affected by optical depth: the to- 
tal gas plus dust density in most of those layers 
(from the disk surface to 1/2 of the disk's height) is 
not high enough to render the disk optically thick. 
This would naturally explain the presence of larger 
grains in the low-temperature component, and is 
consistent with the disk model used for the un- 
derlying continuum which assumes dust grains as 
large as 100 /xm in size in the disk interior (Sect. 
4). 



S 




Fig. 4. — Two dimensional distribution of tem- 
perature of the irradiated accretion disk model 
used to fit the SED of SST-Lup3-1, based on the 
D'Alessio et al. models, illustrating the different 
emission zones. The gas and dust temperatures of 
the disk are given in Kelvin and shown as contours. 
The inset shows in detail the hotter upper layers 
in the inner disk, not visible in the larger figure 
due to their small physical scale. The cold (~70 
K) emission component extends from 5 to 0.7 AU 
and reaches the half disk height in the innermost 
edge. The hot (^300 K) component is mostly lo- 
calized in the inner disk and surface layers up to 
f AU. 

Concerning the crystalline silicates, their mass 
percentages in the hot and cold components are 
16.7% and 12.2% respectively. Given that 59% 
of the dust mass is in the hot component, this 
means that there are ~2 times more crystalline 
grains in the inner disk than in the outer disk. 
This is consistent with the radial distribution of 
crystalline material reported by the high spatial 
resolution ob servations in the inner d isks in Herbig 
Ae/Be stars (jvan Boekel et al.ll2004j ) and suggests 
similar phenomena in disks around the very low 
mass stars. 

Van Boekel et al. 120051 and lApai et all (|2005l ) 



report correlations between high crystalline frac- 
tions and overall large grain sizes or flatter disks 
respectively, based on 10 /jm spectra only. Both 
results suggest a relation between the crystalliza- 
tion and grain growth processes. Our results for 



the hot component alone are consistent with those 
trends, however, the complete 5 to 35 /xm spec- 
trum of SST-Lup3-1 shows a much more complex 
scenario: higher amounts of crystalline silicates in 
the hot component, but not necessarily linked with 
larger amorphous grains, and mostly large amor- 
phous grains in the cold component, where the 
crystalline fraction is only slightly smaller. 

Here we speculate on a scenario implied by the 
observed dust distributions: the inner disk is a 
turbulent region where dust is being partially an- 
nealed when the temperature is above 800 K and 
where grain growth is efficien t for all particles be- 



cause of the higher densities (jKessler-Silacci et al 
2007). The outer disk contains small particles 
with a non-negligible percentage of crystalline ma- 
terial whlchdias either been transported out radi- 
ally (jGaill l2004r ) or which have formed via a dif- 
ferent mechanism than annealing. The 10 /xm fea- 
ture comes mostly from the hot and dense inner 
disk down to a depth in the disk where both large 
and small grains are probed. The long wavelength 
part of the spectrum is emitted by large amor- 
phous grains in deep layers of the disk at distances 
from 0.7 to 5 AU, and by small crystalline grains 
in the disk surface further than 3 AU. According 
to this, the crystalline material could be preferen- 
tially found in the disk surface and in the hot part 
of the innermost disk and the larger amorphous 
grains would be found closer to the midplane once 
grain growth and settling have taken place. Then, 
the high crystalline fractions in disks with large 
grains could be explained if the lower dust opac- 
ities (due to the larger grains) would allow the 
10 /im feature to probe deeper in the inner disk, 
where the l arge and crystalline gra ins are abun- 



dant, as in Ivan Boekel et al. (2005). Such disks 
would also consi stently have a sm aller flaring an- 



gle, as found bv lApai et all (|2005h . On the other 



hand, disks where grain growth has not occurred 
may have small dust with large opacities in the in- 
ner disks, large flaring angles and a 10 (im feature 
coming from the surface of the inner disk, where 
the abundance of large grains is smaller. 

Finally, even assuming the most conservative 
estimate of the disk crystalline mass fraction of 
15%, the disk of SST-Lup3-1 shows a remark- 
ably high degree of crystallinity taking into ac- 
count that only 0.2% of the disk (from 0.01 to 
0.2 AU, out of 10 AU) is warm enough to pro- 
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duce crystalline material. Three main phenom- 
ena may explain this result: i) there is another 
source of disk heatin g such as X-ray flares (e.g., 
Preibisch et al.l [20051 ) or shocks that are operat- 
ing in the system, ii) there was an episodic heat- 
ing event e.g., due to a sudden increase in the 
disk mass accretion rate as obse rved in FU Ori 
stars ( Hartmann fc Kenyonl 19961 ) or iii) there is 
efficient radi al and verti cal transport of dust in 
the system. Gail ( 20041) show that radial mix- 



ing in disks can produce crystallinity fractions of 
several tens of per cent at distances of 5 to 10 
AU, although that model predicts larger amounts 
of forsterite grains in the inner than outer disk, 
contrary to our observations. Future X-ray obser- 
vations of the source, dynamically self-consistent 
disk models of the object, and a larger sample of 
complete 5 to 35 /im IRS spectra of disks may help 
elucidate which is the main effect at work. 

7. Conclusions 

A detailed study of the disk emission of SST- 
Lup3-1 is presented and the results obtained are 
as follows: 

• SST-Lup3-1 is a newly discovered very low 
mass M5.5 young star close to the hydrogen 
burning boundary with a disk, found with 
the Spitzer c2d data in the Lupus III dark 
cloud. 

• The complete 5-35 /im mid-IR Spitzer/IRS 
spectrum of such a low luminosity object is 
presented for the first time. It contains one 
of the highest 10 and 20 /im peak to contin- 
uum ratios reported for such a source. 

• A compositional fit of the spectrum shows 
two temperature components. A hot (T ~ 
300 K) component with a mix of small and 
large grains and a cold component (T ~ 70 
K) with mostly large grains. This latter con- 
clusion depends on the adopted continuum. 

• The spectrum has prominent spectral signa- 
tures indicating significant amounts of crys- 
talline grains in the disk, confirmed in the 
long wavelength range of the spectrum. The 
crystalline mass fraction in the cold compo- 
nent is only slightly less than that in the hot 
disk component. 



• The crystalline grains in the hot component 
are a mix of large (1.5 /im) and small (0.1 
/mi) grains, while those in the cold compo- 
nent are mostly small grains, suggesting that 
the crystalline material is preferentially in 
the disk surface layers. 

• The cold component emission can arise 
from both the "outer" disk surface or from 
deeper layers in the inner disk. In the latter 
case, this would be direct evidence for com- 
bined grain growth and settling with vertical 
height. 



B. M. thanks the Fundacion Ramon Areces for 
financial support. The authors thank Michiel Min 
and to Jeroen Bouwman for fruitful discussions 
on crystalline silicates and the referee for a de- 
tailed report which improved the impact of the 
paper substantially. Support for this work, part 
of the Spitzer Space Telescope Legacy Science 
Program, was provided through Contract Num- 
bers 1256316, 1224608 and 1230780 issued by the 
Jet Propulsion Laboratory, California Institute 
of Technology under NASA contract 1407. As- 
trochemistry at Leiden is supported by a NWO 
Spinoza and NOVA grant, and by the Euro- 
pean Research Training Network "The Origin of 
Planetary Systems" (PLANETS, contract number 
HPRN-CT-2002-00308). 

REFERENCES 

Allard, F., Hauschildt, P. H., & Schweitzer, A. 
2000, ApJ, 539, 366 

Allers, K. N., Kessler-Silacci, J. E., Cieza, L. A., 
& Jaffe, D. T. 2006, ApJ, 644, 364 

Apai, D., Pascucci, I., Bouwman, J., et al. 2005, 
Science, 310, 834 

Apai, D., Pascucci, I., Sterzik, M. F., et al. 2004, 
A&A, 426, L53 

Baraffe, I., Chabrier, C, Allard, F., & Hauschildt, 
P. H. 1998, A&A, 337, 403 

Bouwman, J., Meeus, C, de Koter, A., et al. 2001, 
A&A, 375, 950 

Caballero, J. A., Bejar, V. J. S., Rebolo, R., & 
Zapatero Osorio, M. R. 2004, A&A, 424, 857 



10 



Chapman, N., Evans, II, N. J., Allen, L. E., et al. 
2006, ApJ, submitted 

Comeron, F. 2006, The Lupus Dark Clouds, ed. 

B. Reipurth (Handbook of Low Mass Star For- 
mation in Southern Molecular Clouds), submit- 
ted 

Comeron, F., Fernandez, M., Baraffe, I., 
Neuhauser, R., & Kaas, A. A. 2003, A&A, 406, 
1001 

Comeron, F., Neuhauser, R., & Kaas, A. A. 2000, 
A&A, 359, 269 

D'Alessio, P., Calvet, N., & Hartmann, L. 2001, 
ApJ, 553, 321 

D'Alessio, P., Merin, B., Calvet, N., Hartmann, 
L., & Montesinos, B. 2005, Revista Mexicana 
de Astronomia y Astrofisica, 41, 61 

Dullemond, C. P., Dominik, C, & Natta, A. 2001, 
ApJ, 560, 957 

Evans, II, N. J., Allen, L. E., Blake, G. A., ct al. 
2003, PASP, 115, 965 

Evans, II, N. J., Harvey, P. M., 
Dunham, M. M., et al. 2006, 
http://ssc.spitzer.caltech.edu/legacy/original.html 

Frasca, A., Alcala, J. M., Covino, E., et al. 2003, 
A&A, 405, 149 

Furlan, E., Calvet, N., D'Alessio, P., et al. 2005, 
ApJ, 621, L129 

Gail, H.-P. 2004, A&A, 413, 571 

Guieu, S., Dougados, C, Monin, J.-L., Magnier, 
E., & Martin, E. L. 2006, A&A, 446, 485 

Hartmann, L. & Kenyon, S. J. 1996, ARA&A, 34, 
207 

Honda, M., Kataza, H., Okamoto, Y. K., et al. 
2006, ApJ, 646, 1024 

Jayawardhana, R., Mohanty, S., & Basri, G. 2002, 
ApJ, 578, L141 

Jayawardhana, R., Mohanty, S., & Basri, G. 2003, 
ApJ, 592, 282 

Kessler-Silacci, J., Augereau, J.-C, Dullemond, 

C. P., et al. 2006, ApJ, 639, 275 



Kessler-Silacci, J., Dullemond, C. P., Augereau, 
J.-C, et al. 2007, ApJ, submitted 

Kirkpatrick, J. D., Henry, T. J., & McCarthy, Jr., 
D. W. 1991, ApJS, 77, 417 

Koike, C, Shibai, H., & Tuchiyama, A. 1993, MN- 
RAS, 264, 654 

Lahuis, F., Kessler-Silacci, J., Evans, N., et al. 
2006, c2d Spectroscopy Explanatory Supple- 
ment, Tech. rep., Pasadena: Spitzer Science 
Center 

Luhman, K. L., Briceno, C, Stauffcr, J. R., et al. 
2003, ApJ, 590, 348 

Luhman, K. L., Briceno, C, Rieke, G. H., & Hart- 
mann, L. 1998, ApJ, 493, 909 

Min, M., Hovcnicr, J. W., & dc Koter, A. 2005, 
A&A, 432, 909 

Muzcrollc, J., Hillenbrand, L., Calvet, N., Briceno, 
C, & Hartmann, L. 2003, ApJ, 592, 266 

Natta, A., Testi, L., Muzerolle, J., et al. 2004, 
A&A, 424, 603 

Preibisch, T., McCaughrean, M. J., Grosso, N., 
ct al. 2005, ApJS, 160, 582 

Przygodda, F., van Boekel, R., Abraham, P., ct al. 
2003, A&A, 412, L43 

Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 
618 

Sterzik, M. F., Pascucci, I., Apai, D., van dcr 
Blick, N., & Dullemond, C. P. 2004, A&A, 427, 
245 

van Boekel, R., Min, M., Leinert, C, et al. 2004, 
Nature, 432, 479 

van Boekel, R., Min, M., Waters, L. B. F. M., 
et al. 2005, A&A, 437, 189 

van Boekel, R., Waters, L. B. F. M., Dominik, C, 
et al. 2003, A&A, 400, L21 

White, R. J. & Basri, G. 2003, ApJ, 582, 1109 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



11 



